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TECHNICAL NOTE 58O5 


CALCULATION OP THE FORCES AND MCMENTS ON A 
SLENDER FUSELAGE AND VERTICAL PIN 
PENETRATING LATERAL GUSTS 
^ John M. Eggleston 


SUMMARY 


A theory is presented for calculating the variation with frequency 
of the lateral-force and yawing-moment coefficients due to sinusoidal 
side gusts passing over the profile of a sin5>le fuselage — vertical-fin 
combination. The analysis is based on slender-body theory and is there- 
fore applicable to both subsonic and supersonic airspeeds, provided the 
local flow angles between the profile and the airstream are small. 

The force on any element along the length of a slender body (or flat 
plate) penetrating a gust is found to depend only on the body shape, the 
rate of change of body shape, and the local flow angle. However, when 
the same distribution of local flow angle is obtained by bending of the 
body, the force on 'the body depends both on these same functions and on 
the rate of change of local flow angle along the length of the body. 

The present method of including the penetration effect of the fuse- 
lage and vertical tall in calculating airplane side force and yawing 
moment due to side gusts is believed to be more accurate than the use of 
a simple lag concept to ancount for "the difference in time of penetration 
of the gust by the fuselage and the vertical fin. 


INTRODUCTION 


In calculations involving the lateral forces and moments of an air- 
plane striking a lateral gust, it has been usual to neglect the penetra- 
tion effects and to estimate the forces and moments by considering the 
changes in flow angle to be constant over the configuration. For more 
acctirate analysis of the motion of an airplane passing throu^ gusts, it 
would be desirable to consider the penetration effect and the resulting 
phase variation of the forces and mcments produced by components of the 
airplane which enter the gust at different times. An approximate method 
of accounting for such phase variations due to penetration was described 
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in reference 1, where the effectiveness of several autopilots in contin- 
uous atmospheric turbulence was analyzed. The approximation of reference 1 
involved a linear time lag between the fuselage contribution, Tdiich acted 
at the center of gravity of the airplane, and the tall contribution, which 
acted at the q'uarter-chord point of the mean aerodynamic chord of the tail. 
A more refined method would be to define the pressure or lift distribution 
over the airplane due to a general gust distribution and to confute the 
forces and moments as the gust moves along the body. Such an analysis, 
based on slender -body theory, is presented in this report. 

Slender-body theory was Introduced into the calculation of the forces 
and moments of bodies of revolution and narrow pointed wings in steady 
flow by Munk (ref. 2) and Jones (ref. 3), respectively. Numerous discus- 
sions of the theory, its assunptions, physical concepts, and limitations 
may be found in the literature (see, for instance, the work of Miles in 
ref. it-) and will not be restated herein. 

A generalized solution for a slender body or flat plate in unsteady 
flow was briefly outlined by Garrick in an appendix to reference 5» 3n 
his analysis, Garrick considered a slender airfoil of low aspect ratio 
performing arbitrary vertical translatory or bending motions while passing 
throu^ an Inviscid fluid at uniform velocity. Such an analysis lends 
itself to the Inverse case of a rigid and nonrotating slender lifting 
surface passing at a uniform velocity throu|^ a transverse velocity fluc- 
tuation. using a slender body of revolution and a flat pointed plate 
to represent the fuselage and vertical fin of an airplane, a solution 
for the forces and moments due to the penetration of arbitrary gusts may 
be found, and it is with this particular application that this paper is 
concerned. 

The solution as derived herein is found to be in fairly simple form, 
and the frequency distributions of the lateral forces and moments are 
easily computed. This method of analysis is believed Ijo be more accurate 
than methods which assume no penetration or which assume a linear lag due 
to penetration of an arbitrary gust by an airplane. 

The theory and specific results of this paper may be identified with 
the so-called "Kussner effect." Kussner (ref. 6) derived an expression 
for the lift on a two-dimensional airfoil as it penetrates a sudden 
vertical-gust region without change in direction. The theory was extended 
to finite-span airfoils by Jones in reference 7 "by means of thin-airfoil 
theory. In an appendix to reference 8 are given the Indlclal response 
and sinusoidal-gust response of several wing plan forms of very low aspect 
ratio. Ey properly relating the lifting surfaces analyzed in reference 8 
to represent a fuselage — vertical -fin combination, results for the forces 
due to gusts may be obtained which are in agreement with those obtained 
in the present paper. 
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SYMBOLS 


b wing span 

b pressure altitude 

i 

k reduced -frequency parameter, asc/U 

I lift per unit length dx 

tall length between airplane center of gravity and quarter- 
chord point of mean aerodynamic chord of vertical tall 

p static pressiire 

q dynamic pressure 

s(x) semiwidth of profile along x (see sketch a) 

SqjSj^ dime n sions measured frcm axis of symmetry (see sketch b) 

S wing area of airplane 

t time 

U free-stream velocity 

V transverse or side velocity 

X coordinate in direction of free stream 

Xq,Xi,X 2 fuselage and fin coordinates (see sketch b) 

coordinate normal to free-stream and transverse velocities 
trim angle of attack 
lateral flow angle 
density of air 

airplane relative -density coefficient. 


P 

P 

^^b 


ru 


circular frequency 


Mass 

pSb 



k 
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Mz 

yawing moment about center of gravity 

Py 

side force 

Cn 

yawing-moment coefficient, 

Gy 

side-force coefficient, Py^q,S 

Cl 

lift coefficient based on wing area 

Subscript: 


vt 

vertical tall 


Eie absolute value of a ccsnplex quantity is denoted by j |. 

AHALYSIS 

Lift on a Flat Pointed Plate or Slender Body 
Due to Arbitrary Gustp 

The following sketch shows a flat pointed plate of width 2 s inclined 
to the relative wind at an angle p and subjected to transverse velocity 
v(x,t): 



Sketch (a) 


MCA TN 5805 


5 


In reference 5 ^ I- E. Garrick gives an exjiresslon for the lift per unit 
streamwlse length of a flat plate of increasing area (hut hl^ fineness 
or low aspect ratio) subject to unsteady lift. In terms of v(x,t), the 
expression for the lift per unit length of the flat plate shown in 
sketch (a) becomes: 


Z(x,t) 



z,t)dz 


= irp s 2 (x)~ v(x,t) + aitpU s(x)|| v(x,t) (1) 


where the total derivative, within the approximations of slender-body 
theory, is a function of time and streamwlse distance only: 



+ U 

ox 


( 2 ) 


As derived, then, the lift per unit length is a function of body 
shape s(x), the rate of change of body shape ds/dx, the local transverse 
velocity V, and the rate of change of the local velocity Dv/Dt. Where 
the lifting body is a nonrotatlng flat plate (or of such shape that it 
may be considered a flat plate in two-dimensional transverse flow) and 
the transverse velocity is considered to be some arbitrary gust distri- 
bution, an Interesting slaipliflcation of equation (l) occurs. Since Idie 
gust at any instant of time is a summation of waves mcrving in the 
x-dlrectlon at a uniform velocHy U, the general form of such a gust 
is defined by the relationship 


v(x,t) = v(Ut-x) 


( 3 ) 


It may be obseinred that this relationship is in the form of D'Alembert’s 
classic solution to the general wave equation, which may be found in the 
literature. It is assumed herein that the gust velocity does not vary 
with time as it travels down the lifting body; as in D’Alembert's solu- 
tion, it is assumed that the wave shape does not change with time as it 
moves along the coordinate x. 

For a gust wave, then, 

= 0 
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and ■fche first 1:61711 of equation (l) is likewise zero. Thus, on the flat 
plate or slender body here considered, which is subjected to a gust whose 
velocity does not change with time while it is passing over the body, the 
lift per unit length is a function of the local flow angle only and not 
a function of the local rate of change of flow angle. Ohus 


Z(x,t) = 2itpU s(x)||' v(Ut-x) (ij-) 


Some general ingress ions may be drawn from this development. Since 
the gust velocity does not vary with time as it travels down the lifting 
body, only the shape of the body determines the change of momentum of 
the flow. If the body is an uncambered narrow plate of uniform width or 
a cylinder of imlfoimi diameter, then on these sections of uniform area 
there is no change in momentum of the flow and no forces are developed. 

On the other hand, with a curved cylinder in a steady flow and a strai^t 
cylinder in a tunnel with wavy walls, the momentum of the crossflow would 
change with time and forces would be produced by the lifting surface. 

Where the nose sections of the lifting surfaces just discussed are pointed 
narrow airfoils or slender bodies of revolution, the lift is defined by 
eq'uation (^l- ) . 


Sinusoidal Gust Waves 

Although any harmonic function having the form of equation (3) ml^t 
be employed to represent a gust of arbi-ferary shape and magnitude, probably 
the simplest and most familiar harmonic function is the con^Jlex exponen- 
tial. The gust velocity written in temns of this con5)lex exponential is 
given by 


- . i^(Ut-x) . . 

v(x,t) = U^pe U (5) 


where it may be seen that, for a constant flow angle, v = Up and the 
expressions for lift in the preceding section reduce to the result of 
Jones (ref. 3) for a flat, uncambered airfoil in steady flow. 


Side-Force and Yawing-Moment Coefficients 
for a Simple Airplane Profile 

Since the theory is applicable to airfoils with aspect ratio less 
than 1 and to slender bodies of revolution, a simplified version of the 
general shape of present-day aircraft must be employed. Ihls simplified 
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fuselage — vertical-fin shape is hereinafter referred to as a profile. 

This profile, applicable to a large number of aircraft and also capable 
of treatment by the present theory, is used to compute the side force 
and yawing moment of a fuselage — ^vertical-tail combination. 

nhe profile consists of a nose section of increasing area terminated 
at the center of gravity. That part of the fuselage downstream of the 
center of gravligr will be neglected. Althou^ the theory apply 

to the contracting portion of an isolated fuselage, the flow separation 
effects and the interference due to a wing located in the vicinity of 
the center of gravity make it doubtful that the theory would apply in 
this region. The major contribution to the forces and moments downstream 
of the center of gravl-ty will be considered to be derived frcan the vertical 
fin. Tbus the two components of the simplified profile will be considered 
as Independent lifting sin*faces, any effects due to Interference or cross 
flow being neglected. 

In line with the foregoing discussion, a slnple profile ccanposed of 
a half -ellipse for the forward part of the fuselage and a right triangle 
for the vertical fin is en5>loyed. The major axis of the ellipse is alined 
with the x-axis of the fuselage, and the center of the ellipse is coinci- 
dent with the center of gravity of the airplane. The centroid of the 
rl^t triangle is assumed to be coincident with the center of pressure of 
the vertical fin. 3his profile is shown in the following sketch: 


z 



Sketch (b) 


Qhe shape of the profile is defined analytically by 
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s(x) - Sq 

dfs(x) - sj 
dx 



Xg - 


^X^< X ^ Xg^ (T) 


Exe expression for lift per unit length of the profile Is given by- 
equation (4), and for a sinusoidal latei*al gust defined by equation (5) 
the section lift becomes 


J(x,t) 



2 ii!i(ut-x) 
xe U 


l(x,t) 






K S X S 0) 

< X < *2) 


( 8 ) 


Since the deri-vations of equations (l) and (4) were based on a symmetrical 
plan form, -the contribution of the vertical tail is multiplied by a factor 
of one-half. In doing this, it is assunied that the horizontal tall or 
fuselage acts as an end plate on the half -delta vertical fin, and no other 
correction factor is necessary. 

With the lateral gust referenced to -the center of gravity. 


g(t,x=0) = 3e 


lot 


(9) 


The lift at each streamwise station due to a lateral gust of velocity ^ 
having a unit amplitude is given by the ejqjressions 



-1^ 


(-Xq^ X^ 0) 
(xi< X< Xg) 


> ( 10 ) 


The side-force coefficient and yawing-moment coefficient due to a unit 
side gust are given by the integral over the x dimension of the profile; 
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CyM = ■“ [ -i(x,oj) dx 

^profile 


^[1 -(1 - ii^) 


s 


(11) 


f ^ ^(x,iu) dx 

aoD Jprofile 


r to ~ °o\ 

^2 ” ^ 1 / 


\xoy J__ - - ■ hS\x^ - X. 


= 2£ 
bS 


2 


P^2 . -±^ 

J (x - xi)xe u to 


;:2S^r(2ko - - 12I + fe - =^l)(°l - °n)^ X 

'=0 I- J (k2-kj5 


jjai2 - k^ - 12 - Ik^ka + -^kj^ . I2)e'“^ 


(12) 


■where the notation 


osc 

K, = — -= 
^ U 


(n = 0, 1, 2) 


has been employed to refer to the red,uced frequencies based, on the physical 
dimensions of the profile shown in sketch (b). The mlntis sign appears in 
the expression for the side-force coefficient to agree vith the generally 
accepted conc^t that a positive side gust (moving from starboard to port 
of the airplane) will produce a negative side force at zero frequency. 

Equations (ll) and (l2) may be thou^t of as expressing the in- and 
out-of -phase canponents of the side force and yawing moment due to a unit 
side gust oscillating at some specified frequency. At zero frequency 
equations (ll) and (l2) reduce to ’ 
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CYp = Cy(o>= 0) = - I 2so^ + 


(13) 


^np “ 



m 


NUMERICAL EXAMPLE 


In this section the profiles of three airplanes of various sizes, 
with configurations typical of some current designs, are used to illus- 
trate the variations of side-force and yawing-moment coefficients due 
to unit sinusoidal side gusts of various frequencies. The profile 
dimensions used in the calculations, together with the flight conditions 
and the pertinent airplane dimensions on which the profiles are based, 
are given in table I. 

In applying the theory of this paper to actual aircraft, certain 
practical considerations are involved. Since the ejcpresslons for the 
side force and yawing moment depend only on the configuration, it is 
recognized that wlng-fiiselage interference and other details of a parti- 
cular design may have effects on directional stability which can be 
determined accurately only by wlnd-ttinnel tests or flight measurements. 
Therefore, in order to adjust Cy and to the particular fli^t 
condition chosen, slight modifications in the profiles of the example 
airplanes were made until the theoretical coefficients at zero frequency 
were equal to values of Cy^ and Cup for the example airplanes as 

determined by fll^t tests. This modification consisted in one case of 
a sli^t reduction in tall area and in another case of tripling the fuse- 
lage contribution to adlow for two large engine nacelles and propellers. 
In general, whether these adj\isl3iients are required depends on -Hie air- 
plane configuration. 

The variations with frequency of Cy and due to unit side 

gusts have been plotted in two forms for each of the three exanple air- 
planes; as in-phase and out-of -phase components, and also as the ampli- 
tudes and phase angles relating the unit side gust to the lateral 
coefficients. Hiese variations are shown in figures 1 to 5 for a range 
of frequencies from 0 to 50 radians per second, 

A second set of results is also shown in these figures in order to 
illustrate the difference between the profile theory of this paper and 
the so-called linear i zed-lag theory. This linearized-lag concept, as 
used in reference 1, for example, considers the wing-fuselage contribu- 
tion to act at the center of gravily and the tail contribution to lag 
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TDehind the center of gravity as a function of the time required for the 
airplane to travel a distance equal to the tail length. Expressions 
generally used for this lag are: 


Cyp(o)) « Cyp 


Cnp(o)) « Cnp 


1 - i 


vt 2 ;^! 
CYp IT 




vt 


S “ 


This linear-lag concept has heen applied to the three exarrple airplanes, 
and the results are plotted in figures 1 to 5 in the form, of dashed lines. 

Noted on each figure is the natural frequency of the Dutch roll mode 
of motion of the example aiip)lane. For the particular airplanes consi- 
dered here, the difference in the parameters (side force and yawing 
moment) as obtained by the two methods is small up to approximately twice 
the Dutch roll frequency. At higher frequencies, however, the flow angle 
variation over the airplane has a greater effect, and the linear! zed-lag 
theory is no longer valid. This effect might be important on airplanes 
in -which the Dutch roll frequency is higher in terms of fuselage lengths, 
that is, on airplanes of low density and high directional stability. It 
may be seen -that, for an airplane of short tail length (airplane A, 
fig. 1), both theories are in fair agreement over -the entire range of 
frequencies shown. However, for airplanes of longer -tail lengths, the 
theories agree only at the lower frequencies. Al-though the -varia-tion 
■with frequency of the side force is not generally as Important as -that 
of the yawing moment, it is interesting to note that the side-force 
linearized theory predicts an increasing amplitude with frequency, -whereas 
the present profile theory agrees with physical reasoning that the ampli- 
tude should initially decrease with frequency. 

It is believed, then, that although the frequency -variation of -these 
parameters might differ somewhat -with -the choice of profile shape, the 
trend obtained by considering -the lift distribution over -the profile 
gives results more nearly correct than a theory -which does not consider 
this aspect. In cases -where it is desirable to extend the lateral forces 
and moments to higher frequencies -than the -valid range of the linearlzed- 
lag theory, the present -theory offers a refinement. Since bo-th theories 
require only fairly sin5)le calc-ulations, the compu-tatlon time should not 
be a factor for consideration. 
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CONCLUDING REMARKS 


A simple profile of a fuselage — vertical-fin combination has been 
used to derive a formula for the variation with frequency of the side 
force and yawing moment of an airplane penetrating a gust, and nmerical 
examples of the results are presented. 

In deriving the forces on any element along the length of a slender 
body or flat plate penetrating a gust, it was shown that the force depends 
only on body shape, rate of change of body shape, and the local flcrw 
angle. However, when the same distribution of local flow angle is obtained 
by bending of the body, the force on the body depends both on these same 
fmctions and on the rate of change of local flow angle along the length 
of the body. 

The present method of including the penetration effect of the fuse- 
lage and vertical tall in calculating airplane side force and yawing 
mcmient due to side gusts is believed to be more accurate than the use 
of a simple lag concept to accoxmt for the difference in time of penetra- 
tion of the gust by the fuselage and vertical fin. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., June l4, 1956. 
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TABLE I 

FLIGHT CONDITIONS AND PHYSICAL CHARACTERISTICS 
OP EXAMPLE AIRPLANES 



Airplane A 

Airplane 

Fllgpit Conditions: 

h, ft 

. 50,000 

4,000 

U, ft/sec . 

. 696 

518 

radians 

0.0586 

0.0525 



50 

10.9 



. 0.242 

0.177 

Physical characteristics: 

h, ft 

35.25 

TO 

S, sq ft 

. 250 

540 

S^, sq ft 

55 

71.35 



. 14.8 

27.6 

Xq, ft 

. 10 

15.25 

Xl, ft 

5.7 

22.0 

Sq, ft 

2.7 

3.0 

Sl^ ft 

8.5 

10.6 


Airplane C 
35,000 

700 

0.0825 

31.85 

0.ia^5 

116 

1,1^28 

250 

46.6 

48.5 

56.6 
5.3 

18.5 
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Phase angle, deg 
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Profile theory 

Linearized lag theory 



(Td) Variation with frequency of yawing -moment coefficient 
due to a unit side gust. 


Concluded. 


Phase angle , deg 
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(a) Variation with, frequency of side-force coefficient 
due to a imit side gust. 


Figure 5.- Side-force and yawing -moment coefficients of airplane C as 
detezmined hy profile and linearized theory. 


Phase angle , deg 
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Profile theory 

— — Linearized lag theory 



0 4 8 12 16 20 24 28 

Frequency, o), rodions/sec 


(b) Variation with frequency of yawing-moment coefficient 
due to a unit side gust. 

Figure 5 *- Concluded. 


NACA - Langley Field, V*. 



